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Abstract 



Eigenstates of a particle in a localized and unconfined harmonic potential well are investi- 
gated. Effects due to the variation of the potential parameters as well as certain results from 
asymptotic expansions are discussed. 
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1 Introduction 



Systems confined have received considerable attention in quantum mechanics. The hydrogen 
atom confined in a spherical enclosure was first analyzed in 1937 [I], the restricted rotator 
in 1940 [2] and the harmonic oscillator in 1943 [3]. Since then the one-dimensional oscillator 
immersed in an infinite square potential has received some attention @]-[H]- It should also 
be mentioned that two different cases for a sort of one-dimensional half-oscillator have also 
been reported, one of them bound by an infinite wall [7], [9] and the other one by a finite 
step potential [9] . Recently, the D-dimensional confined harmonic oscillator appeared in the 
literature [10]-[12]. 

The square well potential is an unconfined potential with vertical walls that has been 
used to model band structure in solids [T3] and semiconductor heterostructures [T4]-[T7]. 
The use of a quantum well with sloping sides might be of interest to refine those models. As 
a matter of fact, a sort of localized triangular potential has been an item of recent practical 
[18]-[T9] and theoretical [20]- [23] investigations. 

In this paper we consider the bound-state problem for a particle immersed in an one- 
dimensional harmonic potential which vanishes outside a finite region. To the best of our 
knowledge this sort of trapping has never been solved. By using confluent hypergeometric 
functions the process of solving the Schrodinger equation for the eigenenergies is transmuted 
into the simpler and more efficient process of solving a transcendental equation. Such as for 
the well-known square potential, a graphical method provides some qualitative conclusions 
about the spectrum of this short-range potential well. Approximate analytical results for 
the special cases of low-lying states and high-lying states are obtained with the help of 
asymptotic representations and limiting forms for the confluent hypergeometric function. 
It is shown that the localized harmonic potential yields the full harmonic potential and 
the square well potential as limiting cases. Although the quantization condition has no 
closed form expressions in terms of simpler functions, the exact computation of the allowed 
eigenenergies can be done easily with a root-finding procedure of a symbolic algebra program. 
Proceeding in this way, the whole bound-state spectrum is found. Nevertheless, our purpose 
is to investigate the basic nature of the phenomena without entering into the details involving 
specific applications. In other words, the aim of this paper is to explore a simple system 
which can be of help to see more clearly what is going on into the details of a more specialized 
and complex circumstance such as that one in Ref. [24] . 

2 The particle in a short-range harmonic potential well 

Let us write the short-range harmonic potential well as 

V(x) = Knu 2 (x 2 - rL 2 ) [9 (x + L) - 9 (x - L)} 

{jinu 2 (x 2 — rL 2 ) for \x 
for \x 

where 9 (x) is the Heaviside function, 2L is the range of the potential, m is the mass of the 
particle and u is the classical frequency of the oscillator. The parameter r = \^(0)/[V(0) — 
characterizes four different profiles for the potential as illustrated in Fig. 1. This 
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potential admits scattering states (with E > 0) and bound states (with V(0) < E < and 
r > 0). In what follows we will consider the bound-state problem. 
Let us introduce the new variable 



2mu) 

z = ax, a = \l — ^— (2) 



so that, for |x| < L, the Schrodinger equation 



dx 2 H 2 
turns into the dimensionless form 

d 2 ip(z) ( z 2 



d 2 ib(x) 2m r „ T , . . , , . . 

; + — [E-V(x)]i/j(x) = (3) 



dz 2 \ 4 

where 



+ a)il>(z) = (4) 



V(0) - E 

The general solution for Eq. fll]) can be written as a superposition of definite-parity functions 
[25] 

ip = c e yx(a, z) + c y 2 (a, z) (6) 

where 



/ s , r l a 1 I z 2 \ 



!/ 2 (a, 2 ) = e-^ 2 Af^ + ||0 (7) 

where M(a, 6, 2) = 1-F1 (a; 6; 2;) is the confluent hypergeometric function (Kummer's func- 
tion) 

„ Ml ) = rwf^0f; (8 ) 

1 a 1 (0 + n) n\ 

v ' n=0 v ' 

T (z) is the gamma function, and c e and c Q are arbitrary constants. For x > L, the evanescent 
free-particle solution (ip must vanish as x — )■ 00) is expressed as 

^ce"" 1 (9) 

where c is an arbitrary constant and 



-2mE z 2 ./ a 

Here, 



*=A/nP = Ti/i + ir (10) 



z, = aL=V2^ (11) 
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is the value of z at x = L. Because V (—x) = V (x), the Schrodinger equation is invariant 
under space inversion (x — > —x) and so we can choose solutions with definite parities. The 
even (if) e ) and odd (i[) ) parity eigenf unctions on the entire x-axis can be written as 

Mx) = c e e-^M(^ + ^,^fye(x + L)-9(x-L)] 

+ ce- klxl [6{x-L) + 6{-x-L)] (12) 



ip (x) 



[6(x + L)-0{x- L)\ 



+ ce- fc|x| [9 (x - L) - 9 (-x - L)\ 



(13) 



The even parity solutions satisfy the homogeneous Neumann condition at the origin (dip(x) / dx\ x= o 
0) and the odd ones the homogeneous Dirichlet condition ("0(0) = 0). In this circumstance 
it is enough to concentrate our attention on the positive side of the x-axis and use the con- 
tinuity of ip(x) and dip(x)/dx at x — L. Making use of the recurrence formulas involving y\ 
and y 2 defined in ([7]) [25] 
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one has result 
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The continuity of ^ at x = L says that 
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for even parity solutions, and equal to 
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for odd parity solutions. Matching dip/dx at x — L makes 

- kce- kL 
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for even parity solutions, and equal to 
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for odd parity solutions. Remembering the definition of k from ffTU]) and dividing (1201) by 
(JT7J), and (I2TI) by (ITS]) , one finds the quantization condition 



/ = <? 



(22) 



where 
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for even parity solutions 



for odd parity solutions 
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and 



r a 
4 + 4 



By solving the quantization condition for a in the range 



z L 



< a < 



(24) 



(25) 



one obtains the possible energy levels for a particle trapped in the potential well by inserting 
the allowed values of a in (15]). Hence, 



E = V{0) + \a\ hw 
Notice that a only depends on the potential parameters via r and zl- 



(26) 



3 Qualitative analysis 

A few qualitative results can be obtained with the aid of a plot of the functions / and g on 
the same grid. Figure 2 shows the behaviour of / against \a\ for two different values of y/uL, 
and g for three different values of r. The eigenenergies are determined by the intersections 
of the curves defined by / with the square- root function defined by g (0 < g < y/r/2). 
Without ever solving the quantization condition one is now apt to draw some conclusions 
about the localized oscillator. It is instructive to note that this process for determining 
the spectrum for the localized oscillator looks similar to that one for the square potential. 
Notwithstanding, the zeros and poles of / do not occur at regular intervals as they do for 
tan (x) and cot (x). 
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Seen as a function of \a\ , f presents branches of monotonically increasing curves limited by 
vertical asymptotes due to the zeros of M (a/2 + 1/4, 1/2, z 2 L /2) and M (a/2 + 3/4, 3/2, z 2 L /2). 
For large a/cJL and small \a\, the abscissae of those asymptotes become approximately n+1/2, 
where n is a nonnegative integer, and so do the zeros of /. 

Since the confluent hypergeometric function goes to 1 as zl — > 0, one has that / — > \a\ for 
even parity solutions and / — > — oo for odd ones as y/uJL — > 0. Then, because the square-root 
function vanishes for \a\ = r (zl/2) 2 just one eigenenergy, that one associated with an even 
parity eigenfunction with \a\ ~ 0, is allowed. 

The number of possible bound states grows with r (zl/2) 2 but it is restricted by the value 
of \a\ which makes the square- root function vanish. Therefore, the bound-states solutions 
constitute a finite set of solutions if the potential parameters are finite. 

The spectrum consists of energy levels associated with eigenfunctions of alternate parities. 
The number of allowed bound states grows as the potential parameters increase and there is 
at least one solution, no matter how small the parameters are. All the eigenenergies, in the 
sense of \a\, tend asymptotically to the values n + 1/2 as y/oJL — > oo (n — 0, 1, 2, 3, . . .). The 
energy levels tend to higher energies as the parameter r increases. As a function of \/uL, 
the energy level is a monotonous increasing function for r < 1 but enclosing the oscillator 
with vertical walls (r > 1) makes the energy level to reach a maximum for some value of 

4 Approximate analytical calculations 

Asymptotic representations and limiting forms for the confluent hypergeometric function 
allow us to obtain approximate analytical results for the special cases of low-lying states and 
high-lying states. 

The asymptotic expression which determines M (a, b, z) for a — > — oo reads [25] 

1 _ b 

M{a, b, z) = T{b)e z/2 (hz - azj ' it' 1 ' 2 



X cos 



V2bz - Aaz + Q - - ^ 7T [l + C(|6/2-ar 1/2 )j , for z G R (27) 

Since T (3/2) = T (1/2) /2, it follows that the quantization condition expressed by (122]) takes 
the form 



2 I fv^H-^J t an ^A/jaj^L^ for even parity solutions 
(V\a\z L ) cot (vlajsz) for odd parity solutions 



A further simplification occurs for V(L) — V(0) « \E\ « |^(0)|, when kL » z\/2: 

^J\*A z l^ tan ^A/jof^L^ for even parity solutions 
kL ~ { (29) 
S\a~\z L ^ cot (yV^z^ for odd P arit y solutions 
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In this case the eigenfunction inside the well turns into 



ip (x) ~ c e cos ( y \a\ax I H -= sin ( y \a\ax I (30) 



Here we considered high- lying states in a harmonic potential extending far down (r >> 1) 
and got the solutions for a square well potential. It means that we may neglect any effects 
associated with the bottom of V(x) as far as high- lying states are concerned. It is instructive 
to note that the condition V(L) — V(0) << |^(0)| makes the bottom of the potential look 
fiat. 

For small z, the hypergeometric function M (a, b, z) goes like 

» , , , \ a a (a + 1) , . „ . 

M (aM = l + -z + 1 ^ T ±S + ... (31) 

Thus, the quantization condition for zl « 1 turns into 

I I /, 1 2 r \ 2 n ^ or even parity solutions 
\a\ + [\a\ — j) z L ~ u 



1 _|_ M z 2 ^ Q 

3 L for odd parity solutions 



(32) 

Hence, just one root is allowed: \a\ ~ for the even parity solution. This quasi-null eigenen- 
ergy solution and its very delocalized eigenfunction are valid for y/uL « y^H/m when 
V^(0) ~ 0. In this case the potential looks like a little ripple, a shallow well. 
On the other hand, for large values of \z\ one has [25] 

M(a,b,z) e z z a ~ b 

— — - — ~ — , for Re z > 33 

r (b) r (a) 

In conjunction with the identity T (z + 1) = zT (z) and with the fact that T (z) has simple 
poles at z — —n with n — 0, 1, 2, 3, . . ., the insertion of f )33|) into f l23|) furnishes 

( -1/2 for|a|^n + l/2 
/ * I (34) 
y undefined for | a | = n + 1/2 

for both even and odd parity solutions. The singular behaviour of / when \a\ = n + l/2is the 
reason that it undergoes infinite discontinuities at those values of \a\, as can be grasped from 
Figure 2. It follows that, for sufficiently large Zl, the square-root function can be expressed 
by 



9 * \ (35) 
and the values 

2n + 1/2 for even parity solutions 

(36) 

2n + 3/2 for odd parity solutions 



fulfill the quantization condition. Eq. (1361) represents a convenient approximation as far as 
one considers the lowest values of \a\. As a matter of fact, the intersections of the functions / 
and g occur just slightly below the abscissae of the vertical asymptotes of /. Indeed, a better 
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approximation is obtained as r grows. Nevertheless, for all the values of r, the agreement 
improves as Zl gets larger. In this approximation, M (a, b, z) reduces to a polynomial of 
degree n in z when a = —n. In particular, for 6 = 1/2 and 6 = 3/2 one has [25J 



H 2n 0) 
H 2n +i (x) 



III 



. N „(2n+1)! / 3 9 

n! \ 2 



(37) 



where if n (x) is the Hermite polynomial. Therefore, for y/uiL » y/H/m one gets the 
condensed form 



E n ~ V (0) + 71 + - fc, (X) 



A n e" 



ax \ 



(3f 



where N n is a normalization factor. The approximate results for y/uiL >> y/H/m are 
expected to be exact in the limit L — > oo when the potential goes over to the full-space 
harmonic oscillator. It is comforting to note that the particular values of \a\ obtained from 
the quantization condition are the same as those which make the eigenfunction normalizable 
on the interval (— oo, +oo). The harmonic oscillator approximation for y/uiL finite, though, 
is only reasonable for the low-lying states, i.e. for energy levels so near of the bottom of the 
potential that edge effects can be neglected. 



5 Exact results 

The only remaining question is how to determinate exact results. With the eigenf unctions 
on the whole line expressed by ( fl2l and ( TTBl) . the problem resumes to find the eigenener- 
gies. Although the quantization condition has no closed form solutions in terms of simpler 
functions, the numerical computation of the allowed values of \a\ can be done easily with a 
root-finding procedure of a symbolic algebra program. 

Figure 3 is a plot of the first low-lying energy levels, in the sense of | function 
of y/uiL. The bound-states solutions of the localized and unconfined oscillator constitute a 
finite set of solutions. The number of allowed bound states increases with y/uiL and there is 
at least one solution, no matter how small is y/uiL. All the eigenvalues tend asymptotically 
to the values n + 1/2 as y/uiL — >• oo (n — 0, 1, 2, 3, . . .). The energy levels tend toward higher 
energies as the parameter r increases, as can also be seen in Figure 2. As a function of y/uiL, 
the energy level is a monotonous increasing function for r < 1 but enclosing the oscillator 
with a square well potential (r > 1) makes the energy level reach a maximum for some value 
of y/uiL. 

Figure 4 shows the results for the ground-state eigenfunction against x for y/uiL = 
3/2y/h/m and L equal to a Compton wavelength. Included for comparison is the ground- 
state eigenfunction for the full harmonic oscillator. The normalization J _°° dx |t/>| 2 = 1 was 
done numerically. The eigenf unctions for r = 2 (|a| ~ 0.520) and r = 1/2 (\a\ ~ 0.416) 
differ from that for the full harmonic oscillator. The approximation does better for r = 2, 
as it does for the eigenvalue. In fact, the agreement is not bad even though we have used 
y/uiL ~ y/H/m. Just as expected from the above qualitative analysis, a more successful 
agreement for all the values of r should be obtained for y/uiL >> y/H/m. 
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6 Conclusions 



We have assessed the bound-state solutions of the Schrodinger equation with a localized 
and unconfined harmonic potential well. We have derived the energy eigenvalue equation 
and shown explicitly the eigenfunctions. We have discussed the structure of the solutions of 
the eigenvalue equation. The structure of the eigenfunctions has also been presented. The 
satisfactory completion of this task has been alleviated by the use of graphical methods and 
tabulated properties of the confluent hypergeometric functions. Finally, the exact results 
have been presented. 

As mentioned in the introduction of this work, the three-parameter oscillator potential 
presents a richness of physics which might be relevant for calculations in different fields of 
solid state physics, particularly in electronics and computer components. Furthermore, it 
renders a sharp contrast to the oscillator confined by infinite walls [I]-[E]- 
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Figure 1: Profiles for V(x). 
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Figure 2: Graphical representation for / as a function of \a\. The continuous line for even 
parity solutions, the dashed line for odd ones, and the dotted line for the asymptote. The 
dashed-dotted line stands for the square-root function. 
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Figure 3: | a | for the first three energy levels as a function of ^JZjL (in units of yjh/m) 
for three representative values of r. The dotted parabola stands for the threshold for the 
existence of bound states given by \V(0)\/fku. 
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Figure 4: Eigenfunction for the ground state as a function of x for ^JZjL = and L 

equal to the Compton wavelength. The continuous line for the full harmonic oscillator, the 
dashed line for r = 1/2 and the dotted line for r = 2. 
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